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TIME  AND  SPACE  RESOLVED  ELECTRON  IMPACT  EXCITATION 

RATES 

IN  AN  RF  GLOW  DISCHARGE 

FINAL  REPORT 

A.  Statement  of  Problem  Studied: 

Research  on  rf  glow  discharge  plasmas  was  carried  out  to 
better  understand  the  fundamental  physics  and  chemistry  of 
important  aspects  of  plasma  deposition  and  etching.  A  standard 
reference  system  for  rf  plasma  processing  research  was  utilized. 
Power,  voltage  and  current  waveforms  were  monitored  and  time  and 
space  resolved  measurements  were  made  using  plasma  induced 
emission  (PIE).  Nanosecond  time  resolution  for  PIE  greatly  expanded 
the  utility  of  these  plasma  diagnostics  over  previous  studies. 
Measurements  were  carried  out  with  a  prototype  rf  model  discharge, 
atomic  argon,  concentrating  on  the  important  4s  levels  11.5eV  above 
the  ground  state  and  the  manifold  of  p  levels  coupled  to  these  states 
by  strong  dipole  transitions.  Extensive  analysis  and  computer 
modeling  was  carried  out  to  reproduce  the  experimental  data 
obtained  and  to  test  assumptions  often  used  for  this  model  system. 
Models  of  electron  excitation  "waves"  were  stringently  tested  and 
evaluated.  The  research  concentrated  on  studying  details  of 
excitation  waveforms  as  a  function  of  pressure  and  rf  power. 
Information  on  electron  energy  distribution  functions  were  derived 
from  the  experimental  results.  New  experimental  techniques  to 
enhance  time  and  space  resolution  were  developed. 
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B.  Summary  of  Most  Important  Results: 

In  the  period  1  August  through  31  December,  1991,  significant 
advances  in  the  study  of  electron  kinetics  in  parallel-plate  rf  glow 
discharges  was  made.  The  major  aspect  of  the  investigation  was 
based  on  the  measurement  of  space-  and  time-resolved  optical 
emission  using  time-correlated  single  photon  counting  techniques. 
An  argon  discharge  between  two  equal-area  plane-parallel 
electrodes,  driven  symmetrically  at  13.6  MHz,  was  investigated. 
This  system  closely  approximates  the  ideal  one-dimensional 
discharge  structure  typically  treated  by  modelers  in  the  field.  After 
some  analysis,  experimental  results  can  be  compared  with 
simulation  predictions  of  the  rate  of  electron-impact  excitation. 
This  is  an  essential  test  of  the  treatment  of  electron  kinetics  in  rf 
discharge  models.  In  fact,  the  measurement  of  time-  and  space- 
resolved  emission  spectra  has  been  proposed  as  a  standard 
diagnostic  for  the  GEC  Reference  Cell  as  discussed  at  the  44th 
Annual  GEC,  held  in  Albuquerque,  NM  in  October  1991. 

New  features  of  the  process  of  electron  heating  by  sheath 
oscillation  were  found.1  The  dependence  of  the  range  of  sheath- 
heated  electrons  on  pressure  was  explicitly  revealed  in  space-  and 
time-resolved  emission  data,  (see  Figures  1  and  2)  This  information 
is  important  for  the  verification  of  simulation  results  and  for  the 
interpretation  of  data  acquired  by  other  techniques.  Temporal 
broadening  of  electron-impact  excitation  waveforms  was  also 
measured,  in  agreement  with  simulation  predictions.2  The  effect  of 
finite  spatial  resolution  in  optical  emission  spectroscopy  (OES)  of 
the  parallel-plate  discharge  was  also  illustrated.  Essentially,  the 
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accuracy  of  time-dependent  emission  data  is  limited  by  the  spatial 
resolution  of  the  measurement.  This  problem  results  from  the  3- 
dimensional  nature  of  the  measurement. 

It  is  a  non-trivial  matter  to  determine  the  spatial  resolution 
of  a  typical  imaging  system  used  for  OES.  Normally,  a  portion  of  the 
discharge  is  imaged  onto  the  entrance  slit  of  a  monochromator. 
Previous  investigators  had  found  the  collection  efficiency  of  such  a 
system  for  arbitrary  points  within  a  discharge  volume  by  extensive 
ray  tracing  calculations.4-5  This  method  was  also  used  in  our  initial 
analysis.  However,  the  ray  tracing  analysis  necessary  to  determine 
the  spatial  response  function  requires  a  great  deal  of  programming 
and  computer  time.  We  developed  an  accurate  easy  to  implement 
alternate  method  based  on  the  reasonable  assumption  of  a  nearly- 
ideal  lens  system.6  In  this  approximation,  the  collection  efficiency 
for  an  arbitrary  point  can  be  found  by  geometrical  arguments,  with 
no  need  for  ray  tracing.  A  system  of  achromats,  each  operating  at 
infinite  conjugate  ratio,  easily  satisfies  the  basic  assumptions.  This 
method  reduces  the  time  required  for  programming  and  computation 
by  several  orders  of  magnitude.  The  accuracy  of  this  method  was 
tested  by  comparison  with  exact  ray  tracing  and  the  agreement  was 
excellent.  These  results  should  make  absolute  calibration  and 
spatial  deconvolution  of  emission  data  a  trivial  task  for  researchers 
in  the  field,  (see  details  in  Appendix  1). 

Several  advances  were  also  made  in  the  modeling  aspect  of 
this  project.  A  closer  examination  of  a  collisional-radiative  model 
of  excited  levels  in  an  argon  discharge  led  to  a  number  of  important 
simplifications  which  greatly  expand  the  utility  of  the  model.  The 
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model  was  developed  earlier  as  part  of  a  diagnostic  technique  for 
studying  the  EEDF  in  a  low-pressure  discharge  in  argon.  Using 
absolute  densities  of  the  four  argon  4s  levels,  measured  by  resonant 
laser  absorption,  an  electron  temperature  and  density  consistent 
with  the  model  rate  balance  equations  can  be  found  (a  Maxwellian 
EEDF  is  assumed).  A  simplified  two-level  version  of  this  model  was 
found  to  yield  electron  temperatures  within  10  percent  of  those 
found  using  a  more  complicated  four-level  model.  This  result 
represents  a  great  simplification  in  computation  and  experiment.  For 
the  two-level  version,  the  densities  of  the  lowest  metastable  and 
resonance  levels  must  be  determined.  These  densities  are  typically 
much  higher,  and  therefore,  more  easily  and  accurately  measured, 
than  the  higher  lying  4s  levels.  The  main  advantage  of  this  method  in 
comparison  to  the  emission  line  intensity  ratio  method  is  that  low- 
energy  electrons  are  also  sampled.  An  extension  of  the  model  to 
account  for  a  two-electron-group  EEDF  is  being  considered. 

Progress  was  also  been  made  in  developing  a  "fluid  model" 
code  for  simulating  rf  discharges.  This  type  of  code  is  based  on  the 
first  few  moments  of  the  Boltzmann  equation  for  electrons  and  ions. 
Our  code  takes  into  full  consideration  the  momentum  conservation 
for  electrons  and  ions.  The  solution  is  implemented  by  linearizing 
the  moment  equations  and  adopting  a  simple  implicit  method  of 
discretization.  Preliminary  tests  were  run  on  local  computers  to 
investigate  the  stability  of  the  algorithm  and  to  identify  possible 
areas  of  improvement.  After  further  optimization,  the  code  should  be 
transferred  to  the  facilities  at  the  National  Center  for 
Supercomputing  Applications  (NCSA)  at  the  University  of  Illinois. 
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The  code  will  be  a  powerful  asset  in  guiding  experiments  and 
interpreting  large  amounts  of  experimental  data. 
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APPENDIX 


FAST  CALCULATION  OF  COLLECTION  EFFICIENCY  FOR  OPTICAL 
EMISSION  SPECTROSCOPY  OF  EXTENDED  SOURCES 

M.  J.  Colgan 

Department  of  Physics,  Rutgers  University,  Newark,  New  Jersey  07102 

Abstract 

A  simple  method  for  estimating  the  collection  efficiency  of  a  nearly-ideal  optical  system  is 
described.  Calculation  of  spatial  response  functions  for  the  analysis  of  optical  emission  data 
is  emphasized.  The  simple  geometric  algorithm  presented  is  several  orders  of  magnitude 
faster  than  algorithms  based  on  ray  tracing.  The  method  is  applied  to  the  case  of  a  radially- 
uniform  cylindrical  source  between  parallel  plates  and  the  results  are  in  excellent  agreement 
with  ray  tracing  calculations. 
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FAST  CALCULATION  OF  COLLECTION  EFFICIENCY  FOR  OPTICAL 


EMISSION  SPECTROSCOPY  OF  EXTENDED  SOURCES 

M.  J.  Colgan 

Department  of  Physics,  Rutgers  University,  Newark,  New  Jersey  07102 

I.  INTRODUCTION 

Optical  emission  spectroscopy  (OES)  is  a  standard  technique  for  detecting  atoms  and 
molecules  in  various  excitation  sources  including  glow  discharges,  plasmas,  flames,  arcs, 
etc.  In  many  cases,  the  excitation  source  produces  an  extended  source  of  emission  with 
definite  spatial  structure.  In  order  to  resolve  this  structure,  the  spatial  response  of  the 
measurement  system  must  be  determined  so  that  position-dependent  emission  data  can  be 
deconvoluted. 

A  typical  system  used  for  OES  consists  of  imaging  optics  which  collect  light  from  a 
source  and  focus  the  light  at  the  entrance  slit  of  a  monochromator.  A  photomultiplier  then 
detects  the  light  transmitted  through  the  monochromator.  The  instrument  function  which 
describes  this  system  is  the  product  of  the  efficiency  of  the  detector  (monochromator  plus 
photomultiplier)  and  the  collection  efficiency  of  the  optics  (including  the  entrance  slit  of 
the  monochromator).  If  the  entire  instrument  function  is  known,  position-dependent  emis¬ 
sion  data  can  be  deconvoluted  and  placed  on  an  absolute  scale.  The  detector  efficiency  can 
be  determined  using  standard  calibration  lamps.  The  collection  efficiency  of  the  optical 
system,  however,  can  be  strongly  dependent  on  the  coordinates  of  the  point  of  emission. 
Accurate  calculations  of  collection  efficiency1,2  are  generally  based  on  ray  tracing  algo¬ 
rithms.  While  ray  tracing  is  applicable  to  a  wide  range  of  system  configurations,  the  time 
requirements  for  computer  programming  and  computation  are  extensive.  An  alternative 
to  ray  tracing  is  proposed  in  this  paper  which  is  simple  to  implement  and  reduces  time 
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requirements  by  several  orders  of  magnitude.  The  case  of  a  parallel-plate  discharge  is  con¬ 
sidered  as  an  illustration  of  the  utility  of  these  calculations.  The  results  are  also  compared 
with  ray  tracing  calculations  as  a  test  of  accuracy. 


II.  METHOD 


The  optical  system  under  consideration  is  shown  in  Fig.  1.  The  lens  system  is  arbitrary 
at  this  point.  The  y  axis  coincides  with  the  optical  axis  and  the  z  axis  is  parallel  to  the 
monochromator  entrance  slit.  The  quantity  to  be  determined  is  the  collection  efficiency  e 
of  this  system  for  an  arbitrary  point  (z,  y,  z )  in  the  source  volume.  The  source  is  assumed 
to  be  optically  thin  and  losses  due  to  optical  elements  will  be  neglected.  Isotropic  radiation 
is  also  assumed.  Adopting  the  ray  formalism  for  now,  the  rays  which  are  accepted  by  the 
entrance  slit  will  form  a  pattern  as  they  intersect  a  plane  which  is  tangent  to  the  vertex  of 
the  leftmost  lens  surface  at  y  =  L  (a  convex  surface  is  assumed  for  convenience).  In  most 
cases,  the  ratio  of  L  to  the  effective  aperture  diameter  D  of  the  leftmost  lens  is  greater 
than  about  3.  For  these  conditions,  the  solid  angle  ft  of  light  collection  from  the  point 
(z,  y,  z)  is  very  nearly  equal  to  the  area  A  defined  by  the  pattern  of  accepted  rays  at  y  =  L 
divided  by  (L  —  y)2.  The  collection  efficiency  is  then  given  by 


ft 

£(z,y,z)  =  —  = 


(1) 


4  ir(L  —  y)2 

If  the  lens  system  is  ideal,  i.  e.,  free  of  aberrations,  the  pattern  of  accepted  rays  at 
y  —  L  can  be  found  very  simply.  The  optical  conjugate  of  the  entrance  slit  is  an  accurate, 
possibly  magnified,  two-dimensional  image  of  the  slit  in  the  z-z  plane,  hereafter  referred 
to  as  the  focal  plane.  For  source  points  in  the  focal  plane  and  within  the  confines  of  the 
slit  image,  the  pattern  of  accepted  rays  fills  the  aperture  of  the  leftmost  lens.  For  source 
points  in  the  focal  plane  but  external  to  the  slit  image,  no  rays  are  accepted.  With  these 
assumptions,  the  pattern  of  accepted  rays  from  source  points  which  do  not  lie  in  the  focal 
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plane  can  be  easily  determined  since  the  line  defined  by  the  initial  path  of  an  accepted  ray 
must  pass  through  the  image  of  the  slit.  The  rectangle  defined  by  projecting  those  rays 
which  connect  a  source  point  ( x,y,z )  with  the  comers  of  the  slit  image  onto  the  plane  at 
y  =  L  limits  the  paths  of  rays  which  might  be  accepted.  The  intersection  of  this  rectangle 
with  the  lens  aperture  gives  the  pattern  of  accepted  rays.  The  problem  has  effectively 
been  reduced  to  simple  geometry.  The  parameters  needed  for  solution  are  the  dimensions 
and  location  of  the  slit  image  and  the  diameter  of  the  effective  aperture  of  the  leftmost 
lens.3  This  method  of  determining  collection  efficiency  will  be  referred  to  as  the  projection 
method,  in  contrast  to  the  ray  tracing  method. 

This  ideal  imaging  system  cannot  be  realized  practically.  However,  the  assumptions 
can  be  nearly  satisfied  by  a  system  of  two  computer-optimized  achromats.  Given  the  typ¬ 
ical  dimensions  of  a  spectrometer  entrance  slit,  minimal  spherical  aberration  is  the  most 
important  feature  of  a  nearly-ideal  lens  system.  In  general,  one  of  the  optimization  criteria 
for  a  computer-optimized  achromat  is  that  spherical  aberration  is  reduced  to  the  point 
where  diffraction-limited  focusing  of  parallel  light  is  achieved.  In  this  case,  the  distance 
L  can  be  identified  as  the  “back  focal  length”  of  the  leftmost  lens,  and  is  slightly  shorter 
than  the  focal  length  /.  Note  that  most  achromats  are  asymmetric  elements  and  proper 
orientation  should  be  observed. 

This  degree  of  restriction  on  a  lens  system  may  not  be  acceptable  in  some  situations. 
In  particular,1  cemented  doublets  are  generally  not  useful  for  wavelengths  below  about 
400  nm.  However,  air-spaced  achromats  with  extended  ultraviolet  transmission  ranges  are 
available. 

In  Fig.  2,  a  comparison  of  the  patterns  of  accepted  rays  found  by  the  ray  tracing  and 
projection  methods  is  shown  for  the  source  point  (0.48,0.5,  0.02)  with  coordinates  given  in 
units  of  cm.  For  this  example,  the  lenses  are  identical  /  =  25  cm  achromats  separated  by 
1.0  cm  along  the  optical  axis  and  operating  at  an  //number  of  5.  (The  //number  of  a  lens 


4 


is  defined  as  the  ratio  of  focal  length  to  aperture  diameter.)  The  slit  and  its  image  both 
measure  100  fxm  by  1.0  cm.  Standard  algorithms4  were  used  for  ray  tracing.  The  effects  of 
aberrations  are  apparent  in  Fig.  2,  although  the  areas  enclosed  by  the  two  patterns  differ 
by  only  4  percent. 

The  agreement  between  the  ray  tracing  and  projection  methods  is  generally  very 
good,  with  errors  becoming  significant  only  for  source  points  very  near  the  surface  of  the 
contributing  volume  of  the  discharge.  The  collection  efficiency  for  these  points  is  very  small, 
however,  and  is  difficult  to  calculate  even  by  ray  tracing.  When  the  collection  efficiency 
of  an  optical  system  is  calculated  by  ray  tracing,  the  rays  typically  pass  through  fixed 
grid  points  on  a  lens  aperture  or  some  other  plane.  The  accuracy  for  source  points  with 
small  contributions  is  limited  by  the  spacing  of  this  grid.  This  discrete  nature  of  the  ray 
tracing  method  was  pointed  out  in  Ref.  1.  The  only  way  to  improve  the  accuracy  for  these 
source  points  is  to  decrease  the  grid  spacing  and  trace  a  larger  number  of  rays.  However,  as 
noted  in  Ref.  1,  tracing  even  a  modest  number  of  rays  (8000)  may  require  several  minutes 
of  computer  time  per  source  point.  The  error  due  to  tracing  a  finite  number  of  rays  can 
easily  be  larger  than  the  error  due  to  the  neglect  of  aberrations  in  the  projection  method. 
In  many  applications,  however,  these  sources  of  error  can  be  neglected  based  on  the  low 
collection  efficiency  of  the  source  points  in  question. 

III.  APPLICATION  TO  A  PARALLEL-PLATE  SOURCE 

As  an  example  of  the  utility  of  these  calculations,  the  case  of  a  radially-uniform,  cylin¬ 
drical  parallel-plate  discharge  will  be  considered.  For  this  geometry,  there  are  significant 
variations  in  emission  intensity  along  the  discharge  axis  for  both  DC  and  AC  excitation.5 
Spatially-resolved  emission  spectroscopy  can  determine  excitation,  ionization,  and  dissoci¬ 
ation  rates  yielding  important  information  regarding  local  kinetic  processes  and  discharge 
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sustaining  mechanisms. 

The  discharge  axis  is  assumed  to  coincide  with  the  z  axis  in  Fig.  1.  The  height  of  the 
optical  axis  above  the  lower  electrode  will  be  denoted  by  zo .  The  measured  phown  count 
rate  (single  photon  counting  is  assumed  for  the  detection  mode)  is  then  a(z0)  and  the 
emission  rate  per  unit  volume  is  $(x,y,z  4-  zo)-  The  relation  between  these  two  quantities 
can  be  written  as 


«(*o)  =£o  £(x  ,y,z,z0)$(x,y,z  +  z0)AxAyAz  (2) 

r,y,x 

where  £o  is  the  efficiency  of  the  detector  (including  the  monochromator)  and  e(x,y,z,  z0) 
is  the  collection  efficiency  for  the  source  point  (x,  y,  z  +  zo).  The  sums  over  x,  y,  and 
z  represent  an  integration  over  the  entire  discharge  volume.  The  influence  of  a  nearby 
electrode  on  a  measurement  has  been  included  by  making  £  a  function  of  zq  as  well  as  x, 
y  and  z.  If  radial  variations  are  neglected,  eq.  (2)  can  be  simplified  to 

a(z0)  =  e0^2R(z,z0)$(z  +  z0)Az  (3) 

Z 

where 

tt(z,z0)  =  ^e(x,y,2,*o)AxA*  (4) 

is  the  axial  response  function.  Far  from  either  electrode,  R(z,zo)  is  independent  of  zq  and 
the  notation  R(z,  oo)  will  generally  be  used  unless  the  effect  of  a  nearby  electrode  is  being 
investigated.  As  discussed  earlier,  equation  (3)  represents  the  convolution  of  the  actual 
emission  distribution  and  the  instrument  function  eoR(z,zo). 

This  analysis  would  also  be  applicable  to  the  case  of  time-dependent  emission,  as  in  a 
radio  frequency  parallel- plate  discharge.2,6  In  that  case,  a  and  $  would  be  functions  of  time 
as  well  as  position  in  equations  (2)  and  (3).  The  distortion  caused  by  the  measurement 
system  also  effects  the  time- dependence  of  emission  data  measured  in  these  discharges. 
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In  Ref.  6,  this  effect  was  simply  illustrated  by  reducing  the  //number  of  the  optics.  The 
temporal  widths  of  electron-impact  excitation  waveforms  which  were  derived  from  the 
emission  data  revealed  a  rather  strong  dependence  on  the  //number  of  the  collection 
optics.  The  treatment  considered  here,  which  is  similar  to  that  presented  in  Ref.  2,  is  a 
more  efficient  solution  to  this  problem. 

In  Fig.  3,  R(z,  oo)  is  shown  as  calculated  by  ray  tracing  and  by  the  projection  method 
described  in  the  previous  section.  The  agreement  here  is  excellent.  The  elements  Ax  and  Ay 
were  chosen  so  that  the  variation  in  e  was  small  across  either  distance.  In  this  case,  Ax  was 
set  to  0.02  cm  and  Ay  varied  from  0.02  cm  near  the  discharge  axis  to  0.1  cm  near  the  edges 
of  the  discharge.  R(z,  oo)  was  calculated  for  a  limited  number  of  z  values  by  the  ray  tracing 
method  because  of  the  extensive  computer  time  involved.  For  this  example,  the  discharge 
was  7.62  cm  in  diameter  (x2  +  y2  <  14.2  cm2  in  eq.  (4))  and  the  slit  measured  100  /zm  by 
1.0  cm.  The  lenses  were  both  25  cm  focal  length,  // 5  achromats  (unity  magnification)  and 
were  separated  by  1.0  cm  along  the  optical  axis. 

A  rigorous  series  of  tests  was  performed  to  determine  the  limitations  of  the  projection 
method.  Accuracy  was  not  affected  by  changes  in  the  focal  length  or  //number  of  the 
lenses,  by  magnification  other  than  unity  (using  two  lenses  with  different  focal  lengths),  or 
by  changes  in  the  diameter  of  the  discharge.  Among  the  more  critical  tests  was  the  reduction 
of  the  100  /zm  width  of  the  entrance  slit.  The  agreement  with  ray  tracing  calculations  was 
retained  down  to  a  width  of  20  /zm.  For  a  10  /zm  wide  slit,  a  difference  of  10  percent 
was  observed  at  the  peak  of  the  response  function  ( z  =  0).  However,  the  validity  of  ray 
tracing  for  such  a  narrow  slit  is  questionable  when  diffraction  is  considered.  The  diameter 
of  the  Airy  disk4  for  a  lens  is  2.44  A  f  /D  or  about  6  /zm  for  an  // 5  lens  and  incident 
light  of  wavelength  A  =  500  nm.  The  diffraction-limited  minimum  spot  size  is  thus  nearly 
equal  to  one  of  the  slit  dimensions.  The  projection  method  therefore  has  rougiJy  the  same 
limitations  as  ray  tracing  in  terms  of  a  minimum  slit  width.  Slit  widths  up  to  400  /zm  and 
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slit  lengths  between  0.2  and  2.0  cm  were  also  tested,  demonstrating  agreement  between 
the  two  methods  for  all  cases. 

Optical  systems  are  usually  separated  from  a  source  by  a  window  which  should  be 
considered  as  an  optical  element.  The  primary  effect  of  a  window  is  a  displacement  of  the 
image  along  the  optical  axis.  This  displacement  will  be  very  small  relative  to  the  dimensions 
of  an  extended  source  for  reasonable  glass  thickness.  Quartz  plates  up  to  1.2  cm  thick  were 
added  between  the  discharge  and  the  lens  system  considered  above.  i?(z,  oo)  was  then 
recalculated  by  the  ray  tracing  method  verifying  that  the  plates  had  no  detectable  effect, 
even  when  the  slit  width  was  reduced  to  20  /im. 

Another  parameter  which  might  affect  the  accuracy  of  the  projection  method  is  the 
separation  between  the  two  lenses.  In  this  case,  ray  tracing  showed  a  negligible  change  in 
/2(z,oo)  for  a  separation  of  10  cm  and  a  5  percent  reduction  in  amplitude  for  a  separation 
of  20  cm.  This  small  change  is  probably  due  to  vignetting  and  residual  abberation.  Serious 
errors  might  be  expected  for  very  large  lens  separations. 

The  effect  of  a  nearby  electrode  can  be  investigated  with  a  minor  modification  of  the 
projection  algorithm.  Reflections  from  the  electrode  will  be  neglected  here  as  well  as  the 
curvature  of  the  electrode  edge  facing  the  lens  system.  The  pattern  of  accepted  rays  is 
adjusted  so  that  rays  which  would  strike  the  electrode  are  eliminated  and,  of  course,  there 
is  no  emission  from  points  below  the  electrode  surface.  In  Fig.  4,  the  axial  response  function 
R(z  ,  z0)  is  shown  for  zq  =  oo  and  for  z0  =  0  (optical  axis  in  the  plane  of  the  electrode 
surface).  For  positive  z,  the  ratio  R(z,  0)/R(z,  cc)  is  about  0.5.  This  is  because  the  light 
emitted  from  nearly  all  source  points  with  negative  y  is  blocked  by  the  electrode  when 
zo  =  0.  As  zo  increases,  the  shape  of  R(z,z0)  becomes  rather  complex  before  assuming  the 
shape  of  R(z,  oo).  The  effect  of  any  perfectly  absorbing  obstruction  between  the  source 
and  the  lens  system  can  be  found  in  a  similar  way. 

A  flat,  smooth  reflecting  surface  can  also  be  handled  rather  simply  by  considering  the 
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images  of  the  source  points.  In  the  case  of  a  nearby  plate  electrode  the  image  of  the  source 
point  ( x,y,z  +  z0)  is  (x, y,  —  (z  +  z0 )).  The  collection  efficiency  for  the  image  point  can 
be  calculated  as  if  no  electrode  was  present  and  then  multiplied  by  the  reflectivity  of  the 
electrode  surface.  The  resulting  value  should  then  be  added  to  the  collection  efficiency  for 
the  real  source  point  (x,y,z  +  ^o)- 

Useful  results  can  be  obtained  with  the  projection  method  on  very  short  time  scales. 
While  this  feature  is  generally  appealing,  it  could  be  a  real  advantage  in  situations  where 
the  system  parameters  are  changed  during  the  course  of  an  investigation.  For  example, 
the  spectrometer  slit  width  might  be  reduced  to  improve  the  wavelength  resolution.  As 
noted  earlier,  a  moderate- resolution  ray  tracing  determination  of  collection  efficiency  for 
one  source  point  might  require  several  minutes  of  computation  on  a  personal  computer.1 
In  contrast,  the  results  shown  in  Fig.  3  (solid  line)  were  generated  in  five  minutes  using 
the  projection  method  on  a  personal  computer.  For  that  curve,  collection  efficiencies  for 
more  than  105  source  points  were  calculated. 

IV.  SUMMARY 

A  simple  method  for  estimating  collection  efficiency  and  calculating  spatial  response 
functions  for  optical  emission  spectroscopy  of  extended  sources  has  been  presented.  Re¬ 
sults  for  the  cylindrical  parallel-plate  discharge  geometry  were  compared  with  ray  tracing 
calculations  showing  excellent  agreement  for  a  wide  range  of  system  parameters.  The  con¬ 
ditions  for  which  the  method  is  applicable  place  some  restrictions  on  the  lens  system  which 
may  be  used  for  imaging.  The  fundamental  restriction  is  that  spherical  aberration  must 
be  negligible.  This  condition  is  easily  satisfied  by  a  system  of  two  optimized  achromats, 
each  operating  at  infinite  conjugate  ratio.  Due  to  the  simplicity  and  speed  of  the  algo¬ 
rithm,  the  effects  of  parameter  changes  in  the  optical  system  can  be  rapidly  evaluated. 
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The  response  functions  calculated  in  this  manner  can  be  used  to  remove  the  distortion 
of  position-dependent  emission  data  caused  by  a  measuring  system,  and  to  determine  an 
absolute  scale  for  emission  rates. 
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FIGURE  CAPTIONS 


FIG.  1.  Schematic  of  the  optical  system. 

FIG.  2.  Sample  pattern  of  accepted  rays  at  y  =  L.  The  pattern  of  accepted  rays  determined 
by  the  projection  method  (dashed  line)  is  in  good  agreement  with  the  pattern  found  by 
ray  tracing  (solid  line).  The  areas  enclosed  by  such  patterns  are  used  in  the  calculation 
of  collection  efficiency  and  differ  by  only  4  percent.  These  patterns  were  generated  for  the 
source  point  (0.48,0.5,0.02)  (see  Fig.  1).  For  comparison,  the  pattern  of  accepted  rays  for 
a  point  with  maximum  collection  efficiency,  i.  e.,  the  point  (0,0,0),  would  fill  the  area 
enclosed  by  x2  +  z2  =  2.52. 

FIG.  3.  Axial  response  function  for  a  radially-uniform  cylindrical  source.  The  system  pa¬ 
rameters  are  given  in  the  text.  Nearly  identical  results  are  found  by  the  projection  and  ray 
tracing  methods.  Good  agreement  is  maintained  over  a  wide  range  of  system  parameters. 

Fig.  4.  Effect  of  a  perfectly  absorbing  planar  electrode  on  the  axial  response  function. 
When  the  optical  axis  is  at  the  surface  of  the  electrode  (zo  =0)  there  is  no  source  for 
z  <  0  and  R(z,0)  =  0.  For  z  >  0,  nearly  all  the  emission  from  source  points  with  y  <  0 
is  blocked  by  the  electrode  and  f?(z,  Q)/R(z,  oo)  ~  1/2.  Smooth,  reflecting  planar  surfaces 
can  also  be  treated  simply  using  the  method  described  in  this  paper. 
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